home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / cupgtr.z / cupgtr
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCUUUUPPPPGGGGTTTTRRRR((((3333FFFF))))                                                          CCCCUUUUPPPPGGGGTTTTRRRR((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CUPGTR - generate a complex unitary matrix Q which is defined as the
  10.      product of n-1 elementary reflectors H(i) of order n, as returned by
  11.      CHPTRD using packed storage
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE CUPGTR( UPLO, N, AP, TAU, Q, LDQ, WORK, INFO )
  15.  
  16.          CHARACTER      UPLO
  17.  
  18.          INTEGER        INFO, LDQ, N
  19.  
  20.          COMPLEX        AP( * ), Q( LDQ, * ), TAU( * ), WORK( * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      CUPGTR generates a complex unitary matrix Q which is defined as the
  24.      product of n-1 elementary reflectors H(i) of order n, as returned by
  25.      CHPTRD using packed storage:
  26.  
  27.      if UPLO = 'U', Q = H(n-1) . . . H(2) H(1),
  28.  
  29.      if UPLO = 'L', Q = H(1) H(2) . . . H(n-1).
  30.  
  31.  
  32. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  33.      UPLO    (input) CHARACTER*1
  34.              = 'U': Upper triangular packed storage used in previous call to
  35.              CHPTRD; = 'L': Lower triangular packed storage used in previous
  36.              call to CHPTRD.
  37.  
  38.      N       (input) INTEGER
  39.              The order of the matrix Q. N >= 0.
  40.  
  41.      AP      (input) COMPLEX array, dimension (N*(N+1)/2)
  42.              The vectors which define the elementary reflectors, as returned
  43.              by CHPTRD.
  44.  
  45.      TAU     (input) COMPLEX array, dimension (N-1)
  46.              TAU(i) must contain the scalar factor of the elementary reflector
  47.              H(i), as returned by CHPTRD.
  48.  
  49.      Q       (output) COMPLEX array, dimension (LDQ,N)
  50.              The N-by-N unitary matrix Q.
  51.  
  52.      LDQ     (input) INTEGER
  53.              The leading dimension of the array Q. LDQ >= max(1,N).
  54.  
  55.      WORK    (workspace) COMPLEX array, dimension (N-1)
  56.  
  57.      INFO    (output) INTEGER
  58.              = 0:  successful exit
  59.              < 0:  if INFO = -i, the i-th argument had an illegal value
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.